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Abstract 

A set of generalized squeezed-coherent states for the finite u(2) oscihator is obtained. These 
states are given as hnear combinations of the mode eigenstates with ampUtudes determined by 
matrix elements of exponentials in the su(2) generators. These matrix elements are given in 
the {N + l)-dimensional basis of the finite oscillator eigenstates and are seen to involve 3x3 
matrix multi-orthogonal polynomials Qn(fc) in a discrete variable k which have the Krawtchouk 
and vector-orthogonal polynomials as their building blocks. The algebraic setting allows for 
the characterization of these polynomials and the computation of mean values in the squeezed- 
coherent states. In the limit where N goes to infinity and the discrete oscillator approaches the 
standard harmonic oscillator, the polynomials tend to 2 x 2 matrix orthogonal polynomials and 
the squeezed-coherent states tend to those of the standard oscillator. 

Keywords: Finite quantum oscillator, matrix orthogonality, vector polynomials, Krawtchouk poly- 
nomials, d-orthogonal polynomials, squeezed and coherent states, u(2) algebra. 



1 Introduction 

Discretizations of the standard quantum harmonic oscillator are provided by finite oscillator models 
(see for instance [Ij). We here consider the one based on the Lie algebra u(2) = u(l) © su(2) 
which has been interpreted as a quantum optical system consisting of -|- 1 equally spaced sensor 
points p.J. In this connection, we investigate here the matrix elements of exponentials of linear 
and quadratic expressions in the su(2) generators; these operators represent discrete analogues of 
the squeeze-coherent states operators for the standard quantum oscillator. As shall be seen, these 
matrix elements are given in terms of matrix multi-orthogonal polynomials. 

These polynomials (defined below), generalize the standard orthogonal polynomials by being 
orthogonal with respect to a matrix of functionals [lOj. Very few explicit examples have been 
encountered in the literature; remarkably, our study entails a family of such polynomials and the 
algebraic setting allows for their characterization. 
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1.1 Finite oscillator and u(2) algebra 

The standard one-dimensional quantum oscillator is described by the Heisenberg algebra with 
generators a, and it) obeying 

[a, a^] = id and [a, it)] = [a^ it)] = 0. (1.1) 

The Hamiltonian is given hy H = a'^a + 1/2 and with the position operator Q and the momentum 
operator P defined as follows: 



Q = -(a + at), p = --(a-at), (1.2) 



the equations of motion 



[H,Q] = -iP, (1.3) 
[H,P]=iQ, (1.4) 

are recovered. 

The finite oscillator model is obtained by replacing the Heisenberg algebra by the algebra u(2) = 
u(l) ©5u(2). The su(2) generators are denoted by Ji, J2 and J3 and verify 

[Ji, Jj] = ieijkJk, (1-5) 

with Cijk the Levi-Civita symbol. The u(l) generator is later to be yii). For the finite oscillator, 
the correspondence with the physical "observables" is as follows: 

Position operator: Q = Ji, (1-6) 
Momentum operator: P = —J2, (1-7) 

Hamiltonian: H = J3 + ^ it). (1.8) 

While this relaxes the functional dependence of the Hamiltonian, it is readily seen that this identi- 
fication reproduces the Hamilton-Lie equations (jl.3p and (jl.4p . 

In quantum optics, such a system can be identified with signals coming from an array of + 1 
sensor points |lj. The states of this system can be expanded in the eigenbasis of the Hamiltonian 
H = J3 + N/2 + 1/2, which spans the vector space of the {N + l)-dimensional unitary irreducible 
representation of the su(2) algebra. The eigenstates of H are denoted \N,n) and one has 

H\N,n) = (n+l/2)\N,n), (1.9) 

with n = 0, . . . ,N. The number n will often be referred to as the mode number and the states 
\N,n) as the mode eigenstates. This oscillator model thus only has a finite number of excitations, 
as opposed to an infinite number for the standard oscillator. Moreover, in this representation, the 
spectrum of the momentum and position operators P and Q consists of equally-spaced discrete 
values ranging from —N/2 to N/2. The position and momentum eigenbases can be obtained from 
the mode eigenbasis by simple rotations and their overlaps are su(2) Wigner functions [IJ. 
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It is convenient to introduce the usual shift operators J± and the number operator N . These 
operators are defined by 



J± = (Jl±iJ2), (1.10) 

N = h + N/2. (1.11) 
The action of these operators on the mode eigenstates is given by 



J+\N,n) = yJ{n+l){N -n)\N,n + l), (1.12) 



J-\N,n) = yJn{N -n + l)\N,n-l), (1.13) 
iV|A^,n) = n|iV,n). (1.14) 

For the shift operators J±, the action of any of their positive powers has the form 



/ (n + nMiN — nV , 

Jl\N,n) = ±\N,n + a) = V(-l)"(n + l),(n - N)^ \N, n + a), 

y n\[JS — n — a)\ 



:i.i5) 



J^- \N, n) = ^ ^ - /5) = ^J{-ln-n)^{N-n + l)p \N, n - /3) , (1.16) 

where (n)o = 1 and (n)^ = n(n + 1) • • • (n — a + 1) stands for the Pochhammer symbol. It is 
worth noting that in contradistinction with the standard quantum harmonic oscillator, the finite 
oscillator possesses both a ground state and an anti-ground state. Indeed, one has J+|A^, A^) = 
and J-\N, 0) = 0. This symmetry will play a role in what follows. 

1.2 Contraction to the standard oscillator 

In the limit — t- oo, the finite oscillator tends to the standard quantum harmonic oscillator through 
the contraction of u(2) to t)i [21 [1^. In this limit, after an appropriate rescaling, the shift operators 
J± tend to the operators and a. Precisely, with 

lim iL = a^, lim —= = a, (1.17) 



the commutation relation [a, a^] = id of the Heisenberg-Weyl algebra t)i is recovered. Moreover, the 
contraction of the Hamiltonian H leads to the standard quantum oscillator Hamiltonian //qsc = 
^(P^ + Q^)- This limit shall be used to establish the correspondence with studies associated with 
the standard harmonic oscillator ^12j . 



1.3 Exponential operator and generalized coherent states 

It is known that the standard one-dimensional harmonic oscillator admits the Schrodinger algebra 
5l)i as dynamical algebra [3 This algebra is generated by the linears and bilinears in a and 
, that is a, , it), a^, (a^)^ and a^a. The representation of the group Schi has been recently 
constructed and analyzed in the oscillator state basis by two of us [12j. It involved determining 
the matrix elements of the exponentials of linear and quadratic expressions in a and . The study 
hence had a direct relation to the generalized squeezed-coherent states of the ordinary quantum 
oscillator [9] . 
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We pursue here a similar analysis for the finite oscillator. Notwithstanding the fact that the 
linears and bilinears in J+ and J_ no longer form a Lie algebra, our purpose is to determine 
analogously the matrix elements of the fully disentangled exponential operator 

R{r,, = D{r]) ■ 5(0 = e^'-^+e'^^^e-^'^- • e^-^+Z^e-^-^-Z^, (1.18) 

in the basis of the finite oscillator's states. The parameters r] and ^ are complex-valued and fi = 
log(l + 7]!]). The matrix elements in the (N + l)-dimensional eigenmode basis shall be denoted 

Rk,n = {k,N\R{7j,0\N,n). (1.19) 

In parallel with the definition of the standard harmonic oscillator squeezed-coherent states, we 
introduce the following normalized set of states 

l^'O •= 1/ .1^ ,..^ R{v,O\N,0), (1.20) 
which are a special case of the generalized coherent states 

k 

where ^ is a normalization factor. 

Contrary to the case of the harmonic oscillator, the operator R considered here is not unitary. 
While the operator D^rj) can be shown to be unitary [11], such is not the case for the operator 
S{(,). Nonetheless, we shall observe that the superpositions of states in (|1.20p show spin squeezing 
and entanglement according to the criteria found [6j and [13j. In addition, the consideration of the 
fully disentangled form (ll.lSp allows for the explicit calculation of the matrix elements in terms of 
known polynomials, which is not possible with other choices of the squeezing operator for the finite 
oscillator [15j. 

As previously mentioned, the matrix elements (|1.19p will be naturally expressed in terms of a 
finite family of 3 x 3 multi-orthogonal matrix polynomials Qn{k) in the discrete variable k. In the 
contraction limit, these polynomials tend to the 2x2 matrix orthogonal polynomials encountered 
in pjj. 



1.4 Matrix multi-orthogonality 



Matrix multi-orthogonality has been first studied in the context of Pade-type approximation [5]. 
The algebraic aspects of matrix multi-orthogonality (recurrence relation, Shohat-Favard theorem, 
Darboux transformation, etc.) are discussed by Sorokin and Van Iseghem in |10] . Their study is 
based on matrix orthogonality for vector polynomials. We shall here recall the basic results to be 
used in what follows. 

We first introduce the canonical basis for the vector space of vector polynomials of size q: 













(x\ 



































eo = 




, . . . , e.q—1 — 








, . . . , 625-1 — 




,e2g = 





\xj 



:i.22) 
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For i = Xq + s with i ^ and s = 0, . . . , q — 1, the basis vector has the component in the s + 1*'^ 
position and zeros everywhere else. A vector polynomial of the form aoCo + • • • + an&n with 7^ 
will be said of order n. If g = 1, this corresponds to a standard polynomial of degree n in x. Let 
q and p be positive integers and Hn{x) = (hn,i{x),hn^2{x), . . . ,hn,qix)Y be a (^-vector polynomial 
of order n, where t denotes the transpose [lOj. The vector polynomial Hn{x) is multi-orthogonal if 
there exists a p x q matrix of functionals Qij with i = 0, . . . ,p and j = 0, . . . , g, defined by their 
moments, such that the following relations hold: 



ei,l(/l„,l(x)x'^) + • • • + ei^q{hn,q{x)x'') = 0, 



0, 



,ni - 1, 



(1.23) 



&p,l{hn,l{x)x'') H h Qp,q{hn,q{x)x'') = 0, 



z/ = 0, 



The numbers (ni, 



are defined as follows: set n = fip + (5, with 6 = 0, 



,P 



:i.24) 



1, then 



ni 



n5 



/i + 1 and ns+i 



fl. 



It was shown [10] that such polynomials obey a recurrence relation of the form 

C^^'>Hn+q{x) + ■■■ + C^^^Hn+l{x) + cf 

+ C^-^^Hn-l{x) + • • • + C^-P^Hn-p{x) = xHn{x), 



(1.25) 



along with the initial conditions H^p = ■ ■ ■ = H-i = 0; it was also proven [10\ that a recurrence of 
the type (|1.25p implies the orthogonality conditions (|1.23p and (jl.24p . 

These relations are more easily handled by introducing matrix polynomials, which are matrices 
of polynomials. Suppose that p ^ q, the matrix polynomials are obtained by first writing the 
recurrence relation ()1.25p for k consecutive indices. One has 



x{Hn{x), . . . , Hn+k-l{x)) = {Hn-p{x),... ,Hn+k- 



1+9 



S-p) 



V 



(1.26) 



One can choose k to be the greatest common divisor of p and g; in this case, we can set p = ak, 
q = pk and the matrix on the right hand side can be put in blocks of size k x k. We define the 
qxk matrix polynomial by Qn{x) = {Hnk{x), ■ ■ ■ , Hnk+k-i{x))- The recurrence relation (I1.25P thus 
becomes 



:Q„(x)= rWQ„+,(x). 



(1.27) 



At the end points —a and p in the sum, the matrix coefficient is an upper triangular invertible 
matrix and Tn is a lower triangular invertible matrix. For 5 = 1, this recurrence relation character- 
izes vector orthogonality of order p, also called p-orthogonality. li p = q and F"'^ = (F^)*, matrix 
orthogonality is recovered. In this paper, the special case corresponding to g = 3 and p = 9 will be 
encountered. 
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1.5 Outline 



The outline of the paper is as fohows. In section 2, we obtain the recurrence relation satisfied by 
the matrix elements Rk,n ^-nd show that they involve 3x3 matrix multi-orthogonal polynomials. 
In section 3, we express these matrix elements as a finite convolution involving the Krawtchouk 
polynomials and a family of 3-orthogonal polynomials recently studied in [3]. In section 4, we 
obtain a biorthogonality relation for the matrix polynomials Qn{x). In section 5, we calculate 
the matrix orthogonality functionals Qij for the polynomials. In section 6, we obtain a difference 
equation for the polynomials and discuss the dual picture. In section 7, we derive the generating 
functions and ladder relations. In section 8, we discuss the properties of the states and study 

spin squeezing in this system. In section 9, we briefly review the contraction limit N ^ oo and 
relate our results with those of |12) . We close with concluding remarks in section 10. Appendices 
containing su(2) structure formulas and properties of the Krawtchouk polynomials are included. 



2 Recurrence relation 

We shall begin the analysis by obtaining the recurrence relation satisfied by the matrix elements 
Rk,n- We first observe that 

{k-N/2)Rk,n = {k,N\JsR\N,n) = {k, N\RR-^ JsR\N,n), (2.1) 

where the inverse operator is given by 

= e2^-e-^-^+e^^-e-'''^3g-r?J+_ (2.2) 

The recurrence relation is obtained by expanding the expression R^^J^R and acting on the eigen- 
states \N,n). This is done using the well-known Baker-Campbell-Hausdorff relation 

e^Be-^ = B + [A,B] + ^^ [A, [A, B]] + i [A, [A, [A, B]]] + . . . , (2.3) 

and the formulas for su(2) found in appendix A. Using the polar parametrizations ry = pe''^ and 
^ = re*'*', we obtain 

R-^J^R = (1 - 2p)Ao + pe-'^{l - p)A^ + p(l - p)[e'^ - re-''-^-"'^]A+ 

+ (1 - 2p)re^^^^ - pr^e-'^^~^"'\l-p)Al - 2 pre''^^-^\l - p)A+Ao (2.4) 

Ao = J3+lJl, (2.5) 

= J+-e(l + 2J3)J--f J!. (2.6) 
Introducing this result into (j2.ip . we find that the matrix elements Rk^n obey the recurrence relation 

C^^^Rk,n+3 + C^n^Rk,n+2 + C^^^ Rk,n+1 + C^n^ Rn,k 

+ C^-'^^Rk,n-l + ■■■ + c\^^^Rk,n~9 = k Rk,n- 



where we have defined 

1 + P^ 

A- = J-, 
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(i) 

The coefficients Cn can be obtained straightforwardly with the help of a symbolic computation 
software, their explicit expressions are cumbersome and will thus be omitted here. This recurrence 
relation is of the form (jl.25p with q = 3 and p = 9; consequently, we look for an expression of the 
matrix elements Rk,n as vector polynomials. From the shape of the recurrence relation, it is natural 
to define the 3-vector matrix elements 



^k,n — {Rk,3n, Rk,3n+1, Rk,3n+2) , 

generated by the 3x3 matrix polynomial Qn{k), 

^fc,n = Qn{k)^kfi- 



(2.7) 



(2.8) 



With these definitions, the recurrence relation (|2.7p for the matrix elements Rk^n can be expressed 
as a recurrence relation for the matrix polynomial Qnik)- We have 



kQn{k)= 4^'^ Qn+j{k), 



(2.9) 



where the matrices r„ are expressed in terms of the coefficients Cn in the following manner 



r(i) 



r(-i) 



.(2) 





„(3) 



\ 



-3n+l "-Sn+l 

J-4) J-3) J-2) 

^3n+l ^3n+l ^3n+l 

V-3n+2 '-3n+2 ''3n+2/ 





''3n+2/ 
'^3n > 



r(-2) 



.(1) 

"3x1 



^3n » 



' '-3n 

J-1) JO) 11) 

'^3n+l '^3n+l ''Sn+l 

\''3n+2 ''3n+2 ''3?i+2/ 



^3n 



3n 



'-3n » 
^3n+l 



\''3n+2 '^3)1+2 ^3n+2/ 



(2.10) 



(2.11) 





'-3n 


^3n 


^3n 


r(-3) - 





^(-9) 


^(-8) 




V 





''3n+2/ 



(2.12) 



The structure of the matrix polynomial Qn{k) implies the existence of orthogonality relations of 
the form (|1.23p and (|1.24p . Equivalently, we will explicitly obtain a set of three 3x3 matrices of 
functionals, denoted by J-'i, with orthogonality conditions 



J'i[Qn{k)k'']=0 for i/ = 0,...,L 



n 



(2.13) 



where [xj denotes the integer part of x and where the index i runs from 1 to 3. These orthogonality 
functionals will be computed from the matrix elements of the inverse operator R~^ and will make 
explicit the multi-orthogonal nature of the matrix polynomials Qnik). 
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3 Decomposition of matrix elements 

The matrix elements Rk,n can be expressed as a finite convolution involving the classical Krawt- 
chouk polynomials and a family of vector orthogonal polynomials studied recently in |3J. Indeed, 
one writes 

Rk,n = {k,N\R\N,n), 
N 

= N\ei'^+e^'-^'e-^'^-\N, m){m, N\e^'^+/^e-^'^-^^\N, n), 

m=0 
N 

= ^k,m4'm,n, (3-1) 

m=0 

where we have defined the auxiliary matrix elements 

Xk,m = {k,N\e^^+ef'^'^e-''-^-\N,m), (3.2) 

cpm,n = (m,iV|e€4/2e-?^i/2|Ar,n). (3.3) 

The properties of these intermediary object shall prove useful to further characterize the matrix 
polynomials Qn{k) and will also yield the explicit expansion of the states |?/, ^) in the mode eigen- 
basis. 

3.1 The matrix elements Afc m and Krawtchouk polynomials 

We first study the matrix elements A^^^ of the coherent state operator D{r]) = e'''^+ e^'^^ e~'''^~ . 
These matrix elements can be computed directly by expanding the exponentials in series and using 
the actions (jl.lSp and (I1.16P on the state vector \N,m). With the identity ^^"'^^i = (— n)fc, 
we readily obtain 

m+k i5(k-m) / AJX 1/2 1/2 

2 

where p = and Km{k;p, N) is the Krawtchouk polynomial of degree m, which has the hyper- 
geometric representation 

(— m, — A; 1 \ 
-N > J ' ^^'^^ 

with m = 0, . . . ,N. This result also follows from the recurrence relation satisfied by the matrix 
elements Xkm- This relation can be obtained from (j2.7p or more simply by writing 



(k - N/2) \k,m = {k, N\J3D{7])\N, m) = {k, N\D{7])D~\7])J3D{r])\N, m). (3.6) 
Using the B.-C.-H. relation, the matrix elements Xk,m are found to obey 



+ ^^yJm{N -m + l) \k,m-i- (3.7) 



P 



8 



Introducing the polynomials Xk,m = ^ Km{k;p, N)Xkfi with p = j^^, we recover the 

three-term recurrence relation of the Krawtchouk polynomials 

-kKmik;p,N) = -[p{N -m) + m{l - p)] K^{k;p,N) 

+ p{N - m) K^+i [k- p,N) + m{l- p) K^^i {k; p, N) . (3.8) 

The coherent state operator D{rj) is unitary (see yjj); consequently, we have the following orthog- 
onality relation between the matrix elements Xm,k- 

N 

^k,mXl^n = ^mn, (3.9) 

k=0 

where x* is the complex conjugate of x. In Xk,n, this amounts to the replacement 5 — >• —5. In the 
following, it shall be useful to write the orthogonality relation (I3.9P as the biorthogonality relation 

TV 

y^ Afc,mAjV-n,7V-fc = ^nm, (3.10) 
k=0 

whose equivalence to (|3.9p is shown straightforwardly using the properties of the Krawtchouk poly- 
nomials. Expressing the matrix elements as in (|3.4p yields the well-known orthogonality relation of 
the Krawtchouk polynomials 

(^^^yHl-p)'''''Kr,{k-p,N)K^{k-,p,N) = tlM 5„^. (3.11) 

3.2 The matrix elements (pm,n and vector orthogonal polynomials 

We now turn to the characterization of the matrix elements (pm,n of the squeezing operator S'(^). As 
in the case of the coherent state operator D{r]), the matrix elements (pm.n of S{^) = e^'^+Z^e"^"^-/^ 
can be computed directly by expanding the exponentials in series and applying the actions (ll.lSp 
and (jl.l6p on the state vector \N,n). Obviously, any matrix element (prn'y with m' and n' of 
different parities will be zero; consequently, we set n = 2a + c and m = 2b + c with c = 0, 1 and 
obtain 

6 - ( uaill^rl^^^ l {N-c)\ml / (jV-c)!n! 

'^"''"-^"^^ V {N-m)\ V {N-n)\ ^^'^^^ 

where the identities (0)2™ = 2^" (f)^ (^)n {2(t + s)\ = 2^" s\a\{s + l/2)„ were used and where 
we have defined 

4^)(6;d,Ar) =2^3 (^^5^/2 ^ c^;^)' (3-13) 

with d = — 4r^. The polynomials A^f^ have been studied in [3]; we review here some basic results. 
The polynomials Aa are vector orthogonal polynomials of dimension 3, which corresponds to the 
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case q = 1 and p = 3 of the general setting presented in the introduction. This can be seen by 
computing the recurrence relation satisfied by the matrix elements (j)m,n- Once again, we start with 

{m-N/2)ct>m,n = {m,N\JsS{0\N,n) = {m,N\S{0S-\0J3S{0\N,n). (3.14) 

Using the B.-C.-H. relation and the formulas from Appendix A, we obtain 

S-H^J^SiO = {J3+lJl)+^[J+-ai + 2J3)J--fj^^f. (3.15) 
Substituting this result into (|3.14p yields 



(m - n)4)m,n = {n + l)2(n - A^)2 (f>m,n+2 + (.V (-n)2(A^ - n + 1)2 (/'m,n-2 
3 



with coefficients 



+ ^^E^V^"'')^^'^^ " + (^m,n-2j, (3.16) 

j=0 



/W = {N- 2n)(-l + N + 2Nn - 2n'^), (3.17) 

/(I) = (Gn^ - 12n + N{5 - 6n) + A^^ + 9), (3.18) 

/(2) = (4n-2A^-8), (3.19) 

= 1- (3.20) 



Setting 02b+c,2a+c = ^^-^\J'jfq^;^A^a\b;d,N)(l)2h+c,c, the polynomials A^a\b;d,N) are seen to 
obey the recurrence relation 

(6 - a)A^^{h- d, N) = -^(n + 1)2 (n - N)2Ai%,{b; d, N) - aAi'^b; d, N) 

+ ilY.^-a),fi^^A^:l^{b-d,N). (3.21) 
3=0 

The matrix elements <^m,n are thus given by two families of polynomials Aa for c = 0, 1 of vector 
orthogonal polynomials of order 3. The matrix elements of the inverse operator S~^{£,) can also be 
found by direct computation or by inspection. One readily sees that the matrix elements 0m, n obey 
the biorthogonality relation 

TV 

4'm,n4'*N-n',N-m = ^nn' , (3.22) 

In contrast to the situation with the matrix elements Xk,m of the displacement operator D{i]), the 
biorthogonality relation for the matrix elements 4'm,n is not equivalent to a standard orthogonality 
relation; this is a consequence of the non-unitarity of S{(,) and explains the vector-orthogonal 
nature of the polynomials A^\b;d,N). From the biorthogonality relation of the matrix elements 
(j)m,n follow two biorthogonality relations for the polynomials A^a\b', d, N); we have, for A^ = 2u + 2c, 

|:(-i)'(;)4"(M,~)4i.(«-M,;v) = (i)"*„.. (3^23) 
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For iV = 2^ + 1, we find a biorthogonality relation interlacing the two families c = and c = 1: 

|:(-i)'(;)4'>(M.A'Mf.(»-M,«) = (_„^(.;;„(3^,)^) (i)"*... (3.24) 

3.3 Full matrix elements and squeezed-coherent states 

The results of the two preceding subsections allow to write explicitly the matrix elements Rk,n', 
noting that (j)m,n is automatically zero when m and n have different parities, we have, for n = 2a + c, 
the following expression for the full matrix elements: 

Rk,n = ^ @bK2b+c{k;p,N)A^^\b;d,N), (3.25) 
6=0 

where we have defined 



(^)26+c(^/ ^ xi/2 (jv-c)!(2b + c)! 
°'' = W U + cJ V (^-26-c)! • ^^-^^^ 

The generalized squeezed-coherent states are therefore expressed as the linear combination 

1 ^ 

kOn = ^E^Ml^,fc)> (3-28) 

fe=0 

where Norm is a normalization constant. The expression for the amplitudes simplifies significantly 
if one considers the standard squeezed-coherent states in which the operator R(r],$,) acts on the 
vacuum. Indeed, we have 

If the squeezing parameter ^ = re^"* is set to zero, we recover the standard normalized su(2) coherent 
states 

i''> = \/(rwS(*) 

In section 8, the properties of the states |??,^) will be further investigated; in particular, it will be 
shown that they exhibit spin squeezing. 
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4 Biorthogonality relation 

Given the symmetry of the matrix elements entering the finite convolution yielding Rk,n, the matrix 
elements of the inverse operator S^^{S,)D^^(r]) are expected to have a similar behavior. Indeed, one 
finds that 

N 

Rk,nRN~k,N-n' = ^nn' , (4.1) 

k=0 

where ~ denotes the replacements p — )• —p and r — )• — r. In terms of the vector polynomials 
^k,n = iRk,3n^ Rk,3n+i, Rk,3n+2Y , this biorthogonality relation takes the form 

N 

'ifkA^N-k,N-n'Y = 5„n'Id3x3- (4.2) 

k=0 

This equation can be transformed into a biorthogonality relation for the matrix polynomials Qn{k). 
Indeed, one has 

TV 

Qnik)^k,o{QN-n'{N - k)^N-k,Oy = <5nn'M3x3, (4.3) 

k=0 

which can be written as 

N 

YQnik)W{k){QN-n'{N - k)f = (5wld3x3, (4.4) 
k=0 

with the biorthogonality weight matrix 

(RkflRN-kfl RkflRN-k,l RkfiRN~k,2\ 
Rk,lRN-k,0 Rk,lRN-k,l Rk,lRN-k,2 I • (4-5) 
Rk,2RN-k,0 Rk,2RN-k,l Rk,2RN-k,2/ 

Each element of the weight matrix can be computed exactly as a finite sum with the help of the 
equation (I3.25p . 

5 Matrix Orthogonality functionals 

From the recurrence relation ()2.7p and the results in |10) , it is known that there exists a set of three 
3x3 matrix of functionals with respect to which the matrix polynomials Qn{k) are orthogonal. To 
construct these functionals, we consider the recurrence relation satisfied by the matrix elements of 
the inverse operator i?~^(r7, ^). These matrix elements are defined as 

R-^k = {r^,m~'w,k), (5.1) 

and their recurrence relation is obtained by noting that 

{k-N/2)R-j^ = {n,N\R-^Jo\N,k) = {n,N\R-^JoRR'^\N,k). (5.2) 
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The quantity R ^J^R has been computed previously; recalhng that = J=p and that Jg = J3, the 
recurrence relation for the inverse 3-vector is given by 

3 

kKl= E^n^^^i-,., (5.3) 
i=~i 

where the coefficients are 3x3 matrices and as well respectively lower and 

upper triangular invertible matrices. Since 

N N 

Y,^^]kRKn = Y,{m,N\R-^\N,k){k,N\R\N,n) =5nm, (5.4) 

k=0 k=0 

it follows that 

N 

E'^M(*™!fc)* = Wd3x3. (5.5) 

fc=0 

We may now state the following proposition. 
Proposition 1. The 3-vector can be expressed as 

K!k = ^2P^\k)K.^^^ (5.6) 
i=0 

where Hj = "^^^ . The Ki are 3x3 matrices which depend only on i and p^"'^ (k) are polynomials in 

the variable k. Let n = 3i^ + £, for ^ = 0, 1, 2; if i ^ i, the degree of the polynomial p^[^\k) is v , 
otherwise it is v — \. 

Proof. We set n = in the recurrence relation (|5.3p . which leads to 

= (^M3x3 - e(°))f - e^^-^ - 6(2)^^-1 (5.7) 

(3) 

Since Cn is upper triangular and invertible, it follows that 

%i = ^2^\^^'^^^ll^ (5.8) 

i=0 

where the degree of (A;) is 1 for i = and for all other indices. This establishes ()5.6p for n = 0. 
The proof is then completed by induction. □ 

From this proposition, it is natural to define the following matrices of functionals. 

Definition 5.1. Let J-i for i = 1, 2, 3 be the matrix functionals defined by 

N 

k=0 

With this definition, the relation (|5.5|) can be written as 

n — ? 

Hk'Qnik)] = 03x3 for I. = 0, . . . , , (5.10) 

for i = 1,2, 3. The multi-orthogonality of the matrix polynomials Qn{k) has thus been made explicit 
by the direct construction of the orthogonality functionals. 
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6 Difference equation 



The matrix polynomials Qn{k) are bi-spectral; not only do they satisfy a recurrence relation, they 
also obey a difference equation. In a fashion dual to the approach followed to find the recurrence 
relation, we observe that 

{n - N/2)Rk,n = {k, N\RJ3\N, n) = {k, N\RJ3R-^R\N, n). (6.1) 

Using once again the B.-C.-H relation and formulas from appendix A, we obtain 

RhR-^ = {Bo - ^Bl) - + eH+(l + 2Bo) - ^Blf, (6.2) 

with 

^0 = (1 - 2p)Jo - pe'\l - p)J+ - pe-''{l - p)J_, (6.3) 

B+ = il- p)J+ + 2pe-'^{l - p)J3 - pe-'^'^J_, (6.4) 

H_ = ~pe^'^J+ + 2pe'\l - p)Jo + (1 - p)J-. (6.5) 

Expanding these expressions leads to a difference equation of the form 

6 

nRk,n = X] ^'k^Rk+j,n, (6.6) 

j=-6 

which can be turned into a difference equation for the matrix polynomials Qn{k) which is quite 
involved and not provided here. 

It is interesting to observe that this difference equation contains the same number of terms as 
the recurrence relation, but has a different symmetry. Indeed, the indices run from —6 to 6 in 
the difference equation. It indicates that the order in which the coherent state operator and the 
squeezing operator are presented has an impact on the structure of the associated polynomials. 
If one defines R' = S{^)D{r]) instead of R = D{r])S{^), one gets 6x6 matrix multi-orthogonal 
polynomials satisfying a three-term matrix recurrence relation; these polynomials are not, however, 
orthogonal, because the condition (r^"^-*)* = Fn^ is not fulfilled. 



7 Generating functions and ladder relations 

The ordinary su(2) coherent states \r]) can be used to obtain a generating function for the matrix 
elements Rk,n- We consider the two- variable function defined by 

1 /N\ /N\ 

G(-.rf = (rT7F|^(*) (») 

Clearly, the function G{x,y) can be viewed as the matrix element of R between the coherent states 
|a;) and \y). Introducing a resolution of the identity, we obtain 

N 

Gix,y) = {x\R\y) = ^ (x|I)(r?)|iV, m)(m, iV|5(0|y>. (7.2) 

m=0 
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The first part of this convolution can be evaluated directly using the action of the generators J± on 
the states; not surprisingly, we recover, up to a multiplicative factor, the generating function of the 
Krawtchouk polynomials; the result is 

(x|D(,)|iV, m) = j^^^ (^) "\l + r^xr- {^^m -ij, (7.3) 

2 

with p = -yj^- The second part in the R.H.S of (|7.2p is also readily determined. Setting m = 2t + s, 
one finds 

where z{k) = N — 2k — s. The convolution of the two functions given in (17. 3p and (17. 4p thus 
yields the generating function for the matrix elements The ladder relations for the matrix 

polynomials Qn{k) can also be constructed explicitly from the observations 

V-(n + l)3(n-Ar)3i?fc,„+3 = {k,N\Rjl\N,n), 
^-{-nUN -n + l)^Rk,n-z = {k,N\Rj't\N,n). (7.5) 

The conjugation of the operator by the operator R leads to complicated expressions which are 
best evaluated with the assistance of a computer. 



8 Observables in the squeezed-coherent states 

We now further investigate the properties of the states \r],^) resulting from the application of the 
squeezed-coherent operator D{rj)S{^) on the vacuum |A^, 0). For systems which possess the su(2) 
symmetry, there exist many different parameters to determine whether a state is squeezed or not 
(for a review of the parameters that can be used see [6j). In the following, we will adopt |,13j 

where the indices are to be understood cyclically. If Z^. < 1, the system is squeezed in the direction 
nj, where the nj, i = 1,2,3 are orthogonal unit vectors. This choice of the squeezing criterion is 
relevant because of its relation with entanglement 

The mean value of any observable O in the state can be expressed by 

{O) = 4tE ^^^^^.^y^^' J{2^)\{2j)\{-NU{-N)^, {2j\D-\n)OD{r))\2i), (8.2) 

where k is the normalization constant, which is given by the hypergeometric function 

<r) = mr^hO\'=3Fo (l/2,^,^^;{2r)A . (8.3) 
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For simplicity, we choose to study squeezing along the axis of the Hamiltonian J3. We begin by 

72 ■ 



setting r = to study the behavior of the parameter Z3 in the coherent states. We readily find 



{■h) = iQ) 



— ^cos 5) 
1 + P'^ 



{J2) = -{P) 



Np . 

1 + P'^ 



( J3) = {{H - 1/2)) 



iVl -p2 



{AJ,f = (J|) - {J,f 



Np'' 



a + P 



2\2 ■ 



(8.4) 
(8.5) 



From these results, it is seen that for pure coherent states \r]), we always have Z| = 1, which ensures, 
according to the definition (j8.ip . that purely coherent states are never squeezed; this can be proved 
in a straightforward manner for any choice of normalized basis {ni, 7^2, ^iajllSj. 

We now investigate squeezing in the fi^ = (0,0, 1) direction for the states |?7,0. We have that 



(^1 



pcos((5) 



Gn{r), (J2) 



-psm{5) 

l + p2 



Gn{r) 



{J3) 



l-p 

l + p2 



,21 



Hn{r) (8.6) 



^-^3) = (1^7)2 ([1 + P^] Jn{r) + p^Lnir) + 2p2[cos(25 - 7)]M„(r)) 
where we have defined the polynomials 



^-(^) = 4t E ^(l/2).(-iV)2«(iV - 4i), 



^"(^) = E ^(l/2).(-iV)2.(2z - N/2), 



(8.7) 



(8.9) 



j=0 



LfJ 



Jn{r) = ^ E ^(l/2).(-iV)2.(4i2 - 2iiV + 7VV4), (8.10) 
Kir i] 

1 I 2y 

in(r) = ^ V ^(l/2)i(-iV)2.(-16i2 + 8iiV - iVV2 + N) (8.11) 
Kir i\ 

^ ' 1=0 

= 4t E ^(l/2).+i(-iV)2^+2 (8.12) 
Kir '-^ i\ 

^ ' i=0 

The exact expression for the squeezing parameter cannot be obtained in closed form; nevertheless, 
the squeezing parameter can easily be computed numerically and shows squeezing for various values 
of the parameters; the squeezing parameter is plotted against in figure 1. 
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Figure 1: Squeezing parameter for r = 2, 4, 6 (decreasing amplitudes) with p = 0.8 and = 40. 

9 Contraction to the standard oscillator 

It is of interest to study the behavior of the polynomials Qn{k) and of the squeezed-states in the 
contraction limit N ^ oo where the finite 5u(2) oscillator tends to the standard quantum harmonic 
oscillator. To obtain the proper limit, the parameters of the squeeze-coherent operator must be 
renormalized according to 

O T 

P^-i=, r-^—. (9.1) 

Upon taking the limit, one finds that the 13-term recurrence relation of the Rn,k tends to a 5-term 
symmetric recurrence relation of the form 

2 

kRk,n = Yl Cn^^k,n+j- (9.2) 

i=-2 

This is indeed the type of recurrence relation obtained in [12J. In the case of the standard oscillator, 
it is possible to choose a unitary version of the corresponding S operator and to express its matrix 
elements in terms of 2 x 2 matrix orthogonal polynomials. 



10 Conclusion 

The matrix elements of the exponential operators corresponding to the squeezed-coherent operators 
of the finite oscillator have been determined in the energy eigenbasis of this model. They were 
seen to be given in terms of matrix multi-orthogonal polynomials which have the Krawtchouk and 
vector orthogonal polynomials as building blocks. The algebraic setting allowed to characterize 
these polynomials and to explicitly compute their matrix orthogonality functionals. The results 
have been used to show that the squeezed coherent states of the finite oscillator exhibit squeezing. 
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Appendix A Useful formulas involving the su(2) generators 

The relation 

J3J^ = J^{Js±n), 

holds and can be proven straightforwardly by induction on n. Using this identity and the relations 
(J-t)^ = J=p as well as Jq = J3 , it follows that for P{J±) denoting a polynomial in J±, one has 

[p{J±),Jo] = tJ±p'{J±), 

where P'[x) denotes the derivative with respect to x. The preceding formula and the Baker- 
Campbell-Hausdorff relation lead to the identity 

In addition, we have the relations 

r] = 2nhr-^ + n(n - 1) 

[J_, Jl] = -2nJ^-V3 - n{n - 1) J^'S 

which can also be proved by induction on n. With the help of the previous identities, one obtains 

[J+,P(J_)] = 2JsP'{J-) + J-P"{J-), 
[J-,P{J+)] = -2P'{J+)h - J+P"{J+), 

From these formulas it follows that 

e-P(-^-)j+e-^(-^-) = J+- 2J3P'(J_) - J_[P"(J_) + P'{J-)% 
e^(-^+) J_e-^(-^+) = J_ + 2P'{J+).h + J+[P"{J+) - P'{J+)\ 

A Krawtchouk polynomials 

The Krawtchouk polynomials have the hypergeometric representation 



Kn{x-p,N)=2Fi 



-n, —X 1 
-N 'p 



Their orthogonality relation is 

g (^^^y{l-pr~'^K^{x-,p,N)K^{x-,p,N) = tlM 
They have the generating function 

X N 



nm- 



, =Y.{iyn{x-.p^^r- 

^ ' n=0 ^ ' 



For further details, see [3]. 
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